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report  is  for  accurate  reporting  and  does  not  constitute  an 
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products  or  manufacturers  by  the  National  Aeronautics  and 
Space  Administration. 


Summary 

A  methodology  was  developed  by  Lawrence  W. 
Rehfield  for  the  structural  analysis  of  composite 
rotor  blades.  This  coupled-beam  analysis  is  rela¬ 
tively  simple  to  use  compared  with  the  alternative 
analysis  techniques.  The  beam  analysis  was  devel¬ 
oped  for  thin-wall  single-cell  rotor  structures  and  in¬ 
cludes  the  effects  of  elastic  coupling  achieved  through 
unbalanced  ply  orientation. 

This  paper  demonstrates  the  effectiveness  of  the 
new  composite-beam  analysis  method.  This  is  ac¬ 
complished  through  comparison  of  the  results  of  the 
coupled-beam  analysis  with  those  of  an  established 
baseline  analysis  technique.  The  baseline  analy¬ 
sis  tool  is  an  MSC/NASTRAN  finite-element  model 
built  up  from  anisotropic  shell  elements.  Defor¬ 
mations  for  three  linear  static  load  cases  are  com¬ 
pared.  These  loads  are  centrifugal  force  at  design 
rotor  speed,  applied  torque,  and  lift  for  an  ideal  ro¬ 
tor  in  hover.  A  D-spar  designed  to  twist  under  axial 
loading  is  the  subject  of  the  analysis. 

At  design  rotor  speed  the  finite-element  analy¬ 
sis  and  the  coupled-beam  analysis  indicate,  respec¬ 
tively,  14°  and  15°  of  twist  at  the  spar  tip.  The 
finite-element  model  indicates  less  twist  due  to  rigid 
boundary  conditions  and  wall  thickness  considera¬ 
tions.  A  similar  trend  is  indicated  by  the  applied- 
torque  load  case.  In  the  applied-lift  load  case, 
vertical  deflections  and  twist  indicated  by  both 
analysis  methods  are  essentially  the  same.  Results 
indicate  the  coupled-beam  analysis  is  well  within 
engineering  accuracy. 

The  results  presented  demonstrate  that  moderate 
variations  in  spar  twist  can  be  achieved  by  varying 
the  rotor  rotational  speed.  The  analysis  also  pro¬ 
vides  a  new  and  convenient  approach  for  obtaining 
the  extensional,  torsional,  and  bending  engineering 
stiffnesses. 

Introduction 

A  methodology  was  developed  by  Rehfield  for 
the  structural  analysis  of  composite  rotor  blades 
(ref.  1).  The  beam  force-deformation  relationship  is 
controlled  by  cross-section  properties  that  are  easily 
defined  line  integrals  around  the  rotor  spar.  This 
coupled-beam  analysis  is  relatively  simple  to  use 
compared  with  alternative  analysis  techniques.  The 
beam  analysis  was  developed  for  thin-wall  single¬ 
cell  rotor  structures.  It  includes  the  effects  of 
elastic  coupling  achieved  through  unbalanced  ply 
orientation. 

In  light  of  its  potential,  applications  of  the 
coupled-beam  analysis  are  needed.  One  application 
is  to  design  a  rotor  blade  that  would  change  twist  as 


a  function  of  rotor  speed.  A  scale  model  rotor  blade 
described  in  this  report  was  designed  to  demonstrate 
this  coupling.  The  spar  from  that  preliminary  design 
effort  was  selected  as  the  subject  of  the  analytical 
study  presented  here. 

The  purpose  of  this  paper  is  to  demonstrate  the 
effectiveness  of  the  new  composite-beam  method¬ 
ology.  This  is  accomplished  by  analysis  of  the 
D-spar  through  comparison  of  results  from  use 
of  the  coupled-beam  analysis  and  an  established 
analysis  technique. 

Nomenclature 

By^x  beam  flapwise  curvature,  in”1 

Bz,x  beam  chordwise  curvature,  in"1 

C{j  coupled-beam  stiffness  terms  (from 

ref.  1) 

E\\  lamina  longitudinal  modulus,  psi 

E22  lamina  transverse  modulus,  psi 

EA  beam  uncoupled  extensional  stiff¬ 

ness,  lb 

EIC  beam  uncoupled  chordwise  bending 

stiffness,  lb-in2 

Elf  beam  uncoupled  flapwise  bending 

stiffness,  lb-in2 

G 12  lamina  shear  modulus,  psi 

GK  beam  uncoupled  torsional  stiffness, 

lb-in2 

L  elastic  beam  length,  in. 

Mx  beam  torsional  moment,  in-lb 

My  beam  flapwise  bending  moment, 

in-lb 

Mz  beam  chordwise  bending  moment, 

in-lb 

Nx  beam  axial  force  (centrifugal 

force),  lb 

Qy  beam  shear,  chordwise,  lb 

Qz  beam  shear,  flapwise,  lb 

Sij  coupled-beam  compliance  terms 

(from  ref.  1) 

U,x  beam  extension  strain 

Vx  beam  chordwise  shear  strain 

WjX  beam  flapwise  shear  strain 

w  distributed  beam  load,  lb/in. 

^tip  tip  deflection  due  to  beam  load,  in. 


v\2  lamina  primary  Poisson  ratio 

(j),x  beam  twist  rate,  deg/in. 

Approach 

The  baseline  analysis  tool  selected  is 
MSC/NASTRAN  (refs.  2  to  4).  The  NASTRAN 
finite-element  model  is  built  up  from  anisotropic  shell 
elements.  Deformations  for  three  linear  static  load 
cases  are  compared.  These  loads  are  centrifugal 
force,  ideal  lift,  and  applied  torque. 

The  model  rotor  blade  is  shown  in  figure  1.  The 
spar  from  this  design  (fig.  2)  is  the  subject  of  the 
analysis.  The  finite-element  model  and  the  beam- 
element  model  are  developed  from  the  same  rotor 
geometry.  The  spar  has  a  35.23-in.  radius.  The 
analytical  model  is  rigid  to  radial  station  5.23  in. 
and  has  a  constant  cross  section  from  station  5.23  in. 
to  35.23  in.  The  model  is  divided  into  50  spanwise 
segments  (0.60  in.  each).  The  cross  section  is  defined 
by  16  nodes. 

For  comparison  purposes  the  spar  is  considered 
to  be  a  cantilever  and  is  analyzed  as  a  linear  static 
structure.  It  is  essentially  rigid  from  the  center  of 
rotation  to  the  spar  root.  In  the  finite-element  model 
a  rigid  element  acts  as  an  end  rib  connecting  the 
outboard  end  of  the  beam  element  to  quadrilateral 
elements  at  station  5.23  in. 

The  spar  is  a  high  strain  graphite  and  tough¬ 
ened  epoxy  composite  made  of  Hercules  IM6  fiber 
with  Ciba-Geigy  R6376  resin.  It  has  a  6-ply  layup 
[-1-20/ — 70/-I-20/ — 702/+20],  with  0°  oriented  to  the 
pitch  axis.  The  cured  ply  thickness  is  0.0055  in. 
The  orthotropic  material  properties  are  shown  in 
table  I. 

The  beam  force-deformation  relationship  accord¬ 
ing  to  reference  1  is 


rNxi 

U,x  1 

Qy 

v,x 

Qz 

Mx 

~  \Pij\ 

Wx 

My 

By,X 

Mz\ 

Bz,x\ 

where  [Qy]  is  a  6  by  6  beam  stiffness  matrix.  The 
terms  in  the  stiffness  matrix  are  evaluated  for  the 
D-spar  (table  II).  The  nonzero  coupling  terms  (off- 
diagonal  terms)  are  C14,  C25,  and  C36.  Thus,  cou¬ 
pling  exists  between  extension  and  twist,  between 
flapwise  shear  and  chordwise  bending,  and  between 
chordwise  shear  and  flapwise  bending.  Equation  (1) 
can  then  be  written  as 


NX  =  CUU>X  +  C14^X  lb  (2) 


Qy  =  C22V.X  +  C2^By, x  lb 

(3) 

Qz  =  c33w<x  +  c36bz,x  lb 

(4) 

Mx  =  CUU,X  +  C 44(f) tX  in-lb 

(5) 

My  =  C23VX  +  C33By,x  in-lb 

(6) 

Mz  —  C33W,X  +  CmBz.x  in-lb 

(7) 

In  order  to  apply  forces  and  to  calculate  beam 
deformations,  it  is  necessary  to  invert  equation  (1). 
This  inversion  yields  the  compliance  relationship 


U)X  1 

-Nx- 

y* 

Qy 

Wx 

—  [$ij] 

Qz 

Mx 

By,x 

My 

-  Bz  x  - 

Mz\ 

where  [Qy]  =  [Qy]”1-  For  this  particular  case,  the 
inversion  can  be  accomplished  if  equations  (2)  to  (7) 
are  broken  into  three  pairs  of  equations  with  three 
unknowns  each.  The  following  expressions  for  the 
nonzero  compliance  terms  result: 


Sn  =  b  -  is) 

S22  =  (c*>  -  gj ) 

**  "Kg) 

S44  =  (C44  ~  ^j) 

Ss5  =  (C5S  ■  is) 

5,4  =  -  (is) 511 

S25 = -  (if)  *2 
^ = “  (if) 833 
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Substituting  the  table  II  values  into  equations  (9) 
yields  the  beam  compliance  values  given  in  table  III. 
The  beam  cross-section  force  deformation  can  then 
be  written  in  terms  of  beam  strains  and  curvatures 
as  a  function  of  applied  loads: 


Utx  —  >S"  1 1  Nx  +  SuMx 

(10) 

y,x  =  s^Qy  +  'S25M/ 

(11) 

W,x  =  S33Qz  +  S36Mz 

(12) 

<t>, X  =  SuNx  +  S44MX  deg/in. 

(13) 

By,x  =  5 25Q2/  4“  S33My  in 

(14) 

Bz,x  —  S3qQz  +  SqqMz  in-1 

(15) 

Application  of  Analysis 

Centrifugal  Load  Case 

The  model  rotor  (fig.  1)  has  a  design  rotor  speed 
of  2077  rpm  and  weighs  0.0123  lb/in.  of  span.  These 
values  are  used  to  calculate  the  centrifugal  force  Nx 
applied  to  the  spar: 

Nx  =  934.4  -  0.75287a:2  lb  (16) 

where  x  is  measured  from  station  5.23  in. 

The  material  density  in  the  finite-element  model 
is  selected  to  produce  the  same  weight  per  unit  length 
as  the  model  rotor  (0.0123  lb/in.).  This  is  verified 
by  comparing  the  root  centrifugal  force  calculated 
by  the  finite-element  model  with  the  root  centrifugal 
force  calculated  by  equation  (16). 

Lifting  Load  Case 

The  ideal  rotor  in  hover  has  a  triangular  lift 
distribution  from  radial  station  5.23  in.  to  35.23  in. 
The  total  magnitude  of  the  lift  load  is  selected  to 
be  10  lb.  The  shears  and  moments  produced  by  the 
triangular  lift  and  the  blade  weight  are  calculated  by 
linear  statics.  The  load  w  due  to  lift  and  blade  weight 
can  be  written  in  terms  of  x  as 

w  =  0.02222z  -  0.0123  lb/in.  (17) 

Integrating  once  gives  the  shear 

Qz  =  0.01111a:2  -  0.0123a:  -  9.631  lb  (18) 
Integrating  again  gives  the  moment 

My  =  0.003704a:3  -  0.00615a;2 

-  9.631a;  +  194.5  in-lb  (19) 


In  the  case  of  the  finite-element  model,  the  lift 
is  applied  as  concentrated  forces.  The  forces  are 
applied  vertically  at  the  spanwise  row  of  grid  points 
along  the  upper  surface  and  at  the  quarter-chord. 

Torsional  Load  Case 

The  torsional  load  (torque)  is  applied  at  the  tip 
of  the  spar.  The  torque  is  constant  with  respect  to 
length  and  is  selected  to  be 

Mx  =  50.0  in-lb  (20) 

For  the  finite-element  model,  the  torque  is  divided 
equally  among  the  nodes  at  the  tip  and  is  applied 
as  concentrated  moments.  An  end  rib  composed  of 
shear  elements  is  used  in  this  load  case  to  prevent 
excessive  distortion  of  the  cross  section.  These  ele¬ 
ments  provide  no  resistance  to  warping  in  the  axial 
direction. 

Calculation  of  Deflections  for  Coupled-Beam 

Analysis 

The  internal  beam  forces  (centrifugal  forces, 
shears,  and  moments)  are  calculated  for  the  spanwise 
stations  by  means  of  equations  (16)  to  (20)  and  av¬ 
eraged  over  the  segment  length.  The  forces  are  then 
multiplied  by  the  beam  compliance  matrix  to  yield 
beam  deformation  derivatives.  These  derivatives  are 
integrated  along  the  span  to  produce  deflections. 

Results 

Centrifugal  Load  Case 

In  the  centrifugal  load  case,  the  coupled-beam- 
element  analysis  and  the  finite-element  model  agree 
(fig.  3).  The  coupled-beam  analysis  predicts  —15.29° 
of  twist  at  the  tip,  and  the  finite-element  model 
predicts  -14.13°  at  the  tip.  Figure  4  depicts  plotter 
output  from  the  finite-element  analysis  showing  this 
twist. 

The  difference  in  the  twist  indicated  by  the  two 
analyses  can  be  attributed  to  several  effects.  The 
boundary  effect  induced  by  the  rigid  inboard  end  rib 
has  a  stiffening  effect  on  the  finite-element  model. 
The  rigid  element  does  not  permit  any  warping  of 
the  cross  section  out  of  its  plane.  This  locally  stiffens 
the  structure  in  torsion.  The  coupled-beam-element 
model  used  in  this  case  assumes  free  warping  at  the 
root.  Additionally,  the  rigid  element  does  not  permit 
Poisson  contraction  at  the  root.  This  has  a  local 
stiffening  effect  in  extension,  and  thus  reduces  the 
twist  associated  with  extension.  The  coupled-beam 
analysis  assumes  free  Poisson  contraction.  These 
local  root  stiffening  effects  are  apparent  in  figure  5. 
The  effects  begin  to  die  out  at  approximately  radial 
station  10  in. 
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The  coupled  response  of  the  spar,  outboard  of 
radial  station  10  in.,  shows  good  agreement  between 
the  two  analyses.  (See  fig.  5.)  The  constant  offset 
between  the  two  curves  can  most  likely  be  attributed 
to  wall  thickness  considerations.  The  coupled-beam 
analysis  is  a  thin-wall  analysis  that  does  not  include 
the  composite-plate  torsional  stiffness.  The  finite- 
element  model  includes  this  stiffness. 

Lifting  Load  Case 

The  two  analyses  also  agree  well  on  vertical  de¬ 
flections  due  to  combined  lift  and  blade  weight.  The 
coupled-beam  analysis  predicts  6.70  in.  of  deflection 
at  the  spar  tip.  The  finite-element  analysis  predicts 
6.48  in.  of  deflection  at  the  tip.  (See  fig.  6.)  The 
deflection  predicted  by  the  coupled-beam  analysis 
is  greater  than  that  predicted  by  the  finite-element 
analysis. 

The  coupled-beam  analysis  indicates  that  spar 
twist  is  uncoupled  from  both  shear  and  bending 
(fig.  7).  The  finite-element  analysis  indicates  a  small 
coupling  that  can  be  considered  as  insignificant  for 
design  purposes  (fig.  7). 

Torsional  Load  Case 

In  the  applied-torque  case  agreement  is  also  good. 
Through  use  of  equations  (13)  and  (20),  twist  <j)  at 
the  spar  tip  can  easily  be  expressed  for  Nx  =  0  as 
follows: 

(j)  =  S44MxL(180/tt) 

=  (0.2003  X  10-3)  (50.0)  (30.0)  (180/tt) 

=  17.2° 

The  twist  predicted  by  the  finite-element  model 
is  distorted  at  the  tip  because  of  the  concentrated 
moments  at  the  end  nodes.  (See  fig.  8.)  The  apparent 
excellent  agreement  of  the  two  analyses  at  the  tip  is 
due  to  this  distortion.  The  true  agreement  between 
the  methods  is  better  shown  in  the  twist  rate  curve. 
(See  fig.  9.)  Again,  the  constant  offset  between  the 
two  curves  in  figure  9  can  most  likely  be  attributed 
to  wall  thickness  considerations. 

Comparison  to  Engineering  Beam  Theory 

Bernoulli  beam  theory  gives  us  the  uncoupled 
relationship  between  beam  curvature  J3X,  bending 
moment  M,  and  bending  stiffness  El  as 

B,x  =  M/ El 

A  similar  statement  can  be  made  for  the  compos¬ 
ite  beam  if  coupling  effects  are  ignored.  If  the 


beam  chordwise  shear  strain  is  assumed  to  be  zero, 
equation  (6)  becomes 


Then  C55  corresponds  to  the  beam  uncoupled  flap- 
wise  bending  stiffness  Elf.  Alternately,  if  the  chord- 
wise  shear  force  is  taken  to  be  zero,  equation  (14) 
becomes 

By,x  —  S55  My 

In  this  instance,  the  inverse  of  S55  corresponds  to 
EIj.  In  the  case  of  the  cantilever  beam,  beam 
strains  are  not  constrained  and  bending  deflections 
due  to  applied  beam  loads  are  sought.  For  this 
case  it  is  correct  to  take  the  beam  stiffness  from  the 
compliance  term  (S55)  rather  than  from  the  stiffness 
matrix  term  (C55).  The  significance  of  selecting  the 
correct  term  is  illustrated  by  the  particular  case. 

For  a  uniform  beam  with  a  triangular  load,  the 
deflection  at  the  tip  (ref.  5)  is  given  as 

ztip  =  11  wtipL4/l20EIf  in. 

where  u>tjp  is  the  maximum  height  of  the  triangular 
lift  distribution.  For  the  beam  weight  alone,  the  tip 
deflection  is  (ref.  5) 

zt[ p  =  —wL^/SEIf  in. 

Superimposing  the  two  cases  and  considering  the 
10-lb  triangular  lift  load  and  given  blade  weight  of 
0.0123  lb/in.  results  in 

zt[p  =  48  254.6 /Elf  in.  (21) 

In  our  particular  case,  C55  is  almost  twice  as  large 
as  St^1 .  Since  it  is  not  a  realistic  constraint  to  force 
the  chordwise  shear  strain  to  zero,  the  beam  stiffness 
should  not  be  taken  to  be  C55.  For  our  case  the 
chordwise  shear  force  is  zero,  so  Elf  is  taken  to  be 

S"1.  Substituting  this  value  into  equation  (21)  yields 
zt\p  =  6.544  in. 

This  value  agrees  well  with  the  deflection  predicted 
by  Bernoulli  beam  theory  and  by  finite-element  anal¬ 
ysis.  The  difference  is  primarily  due  to  rounding 
off  of  the  constant  in  the  tabulated  solution  of  ref¬ 
erence  5.  Taking  the  flapwise  bending  stiffness  to 
be  C55  would  incorrectly  indicate  that  the  beam  is 
twice  as  stiff.  Similar  reasoning  shows  that  since  the 
flapwise  shear  strain  is  not  constrained,  the  beam 
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chordwise  bending  stiffness  EIC  is  taken  from  the 
compliance  relationship  as  ■W  • 

Since  beam  extension  is  not  constrained,  the  tor¬ 
sional  stiffness  GK  is  taken  from  the  compliance  re¬ 
lationship  as  ^44  ‘  Similarly,  since  beam  twist  is  not 
constrained,  the  extensional  stiffness  EA  is  . 

The  engineering  stiffness  constants  are  summa¬ 
rized  in  table  IV.  It  should  be  emphasized  that  these 
stiffnesses  are,  in  general,  dependent  on  boundary 
conditions  and  are  therefore  not  simply  cross-section 
properties. 

Conclusions 

Results  from  the  coupled-beam  theory  agree 
well  with  the  finite-element  analysis  predictions. 
At  design  rotor  speed  the  finite-element  analysis 
and  the  coupled-beam  analysis  indicate,  respec¬ 
tively,  14°  and  15°  of  twist  at  the  spar  tip.  The 
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Table  I.  Orthotropic  Material  Properties 


-Eli,  psi 
#22,psi 
^12  • 
Gi2,psi 


23.1  x  106 
1.4  x  106 
.  .  0.338 
0.73  x  106 


Table  II.  Nonzero  Beam  Stiffness  Terms 


On, lb  .  . 
0*22,  lb  • 
C33,  lb  . 
C44,  lb-in2 
O55,  lb-in2 
066,  lb-in2 
O14,  lb-in. 
O25,  lb-in. 
036,  lb-in. 


0.8332  x  106 
0.1651  x  106 
0.3071  X  105 
.  0.9747  x  104 
.  0.1337  x  105 
.  0.1128  x  106 
-0.6294  X  105 
0.3147  x  105 
-0.3147  X  105 


Table  III.  Nonzero  Beam  Compliance  Terms 


-Sn,  lb"1  .  . 
S22,  lb”1  •  • 
533,  lb”1  •  • 

S44,  (lb-in2) 
555,  (lb- in2) 

56e,  (lb-in2) 
S14,  (lb-in.)"1 
525,  (lb-in.)-1 

536,  (lb-in.)-1 


0.2345  X  10"5 
0.1099  x  10-4 
0.4561  x  10-4 

0.2003  x  10"3 
0.1356  x  lO-3 

0.1242  x  10-4 

0.1513  x  10-4 
-  0.2585  x  10-4 
0.1273  x  10-4 
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Table  IV.  Summary  of  Beam  Uncoupled  Engineering  Stiffness  Constants 


7.5  r 


-  Coupled-beam  analysis 

- Finite-element  analysis 


figure  4.  Finite-element  model  twist  due  to  centrifugal  force 
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Figure  5.  Twist  rate  due  to  centrifugal  force. 


Figure  6.  Beam  deflection  due  to  lift  and  blade  weight. 
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Figure  8.  Twist  due  to  applied  torque. 
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